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摘 要 : 利用 延 河 流域 1998 年 .2004 年 .2010 年 和 2016 年 4 虹 


的 Landsat 遥感 影像 ,计算 遥感 生态 指数 (Remote 


Sensing Ecological Index, RSET) ,对 流域 生态 环境 质量 的 时 空 变化 进行 分 析 和 评价 ,利用 CA-Markov 模 型 对 2022 年 
和 2028 年 的 生态 环境 进行 模拟 和 预测 。 结 果 表 明 :(1) RSEI 评 价 指数 在 延 河流 域 具 有 较 好 的 适用 性 ,可 用 于 监测 
和 评价 流域 生态 环境 的 时 空 变化 特征 。(2) 时 间 尺 度 上 ,RSFEI 均 值 从 1998 年 的 0.392 变化 为 2016 年 的 0.530, 呈 逐年 
上 升 趋势 ,空间 尺度 上 ,RSEI 从 东南 向 西北 逐渐 增 大 。 以 0.2 为 间隔 划分 RSEI 等 级 ,发现 流 域 生态 环境 质量 等 级 以 
较 差 和 中 等 为 主 , 占 总 面积 的 62.4% ,主要 分 布 于 流域 的 山谷 地 带 , 优 良 等 级 的 面积 占 比 最 小 , 约 为 8.1% ,主要 分 布 
于 流域 南部 的 山区 地 带 。(3) 自然 和 人 为 因素 对 流域 生态 环境 质量 变化 的 影响 存在 较 大 空间 差异 性 。 其 中 海拔 在 


1200~1800m 的 区 域 , 生 态 环境 质量 变化 主要 受 自 然 因 素 影 响 ; 而 海拔 低 于 1200m 的 区 域 , 则 主要 受 人 类 活动 影 
响 。(4) 预测 到 2028 年 ,流域 生态 环境 质量 会 持续 提升 ,良好 和 优良 等 级 的 面积 占 比 将 达到 66% ,但 是 流域 西北 部 


由 于 生态 环境 本 底 脆 弱 , 仍 存在 较 大 的 问题 。 


关键 词 : 延 河 流域 ;， RSEI; 生态 环境 质量 ; 主 成 分 分 析 ; CA-Markov 模 型 


生态 环境 作为 人 类 生存 的 基本 保障 和 社会 赖 
以 发 展 的 物质 基础 ,其 状态 会 直接 影响 人 类 生活 的 
质量 "。 党 的 十 九 大 报告 将 “美丽 ”作为 社会 主义 现 
代 化 强国 的 基本 特征 ,把 “坚持 人 与 自然 和 谐 共 生 ” 
纳入 新 时 代 坚 持 和 发 展 中 国 特 色 社 会 主义 方略 。 
在 未 来 ,我 国 将 继续 重点 实行 经 济 与 生态 环境 可 持 
续 发 展 战略 ,大 力 加 强生 态 环境 监测 工作 ”。 

生态 环境 质量 测评 方法 主要 包括 两 类 :一 类 是 
单 因子 变化 分 析 ,例如 分 析 土 地 利用 变化 ”植被 净 
初级 生产 力 变化 “和 植被 覆盖 度 变化 ”等 与 生态 环 
境 密 切 相 关 的 因子 变化 。 男 一 类 是 多 因子 变化 综 
合 分 析 , 与 单 因子 分 析 相 比 ,多 因子 综合 分 析 更 加 
全 面 和 准确 ,学 者 们 也 为 此 提出 了 多 种 评价 指标 体 
系 “" 。 有 目前 ,应 用 最 广泛 的 是 我 国 于 2006 年 发 布 
的 《生态 环境 状况 评价 技术 规范 》 中 提出 的 EI(Eco- 
logical Index ) 指数 ”和 徐 涵 秋 ”于 2013 年 提出 的 
RSEI(Remote Sensing Ecological Index ) 指 数 。 

RSEI 是 将 可 以 反映 生态 问题 的 绿 度 指 标 (ND- 


VI) 湿度 指标 (Wet) .十 度 指标 (NDBSI) 和 热度 指标 
(LST) 通 过 主 成 分 分 析 方 法 进行 耦合 ,是 一 个 可 以 量 
化 的 生态 质量 指数 。 相 比 于 EI, 它 解决 了 权重 设置 
合理 性 问题 和 生态 环境 质量 可 视 化 问题 "。 目 前 ， 
RSEI 多 用 于 评价 城市 和 矿区 的 生态 环境 质量 "7 ， 
而 针对 流域 系统 开展 生态 环境 质量 评价 研究 相对 
较 少 。 同 时 ,对 影响 生态 环境 质量 的 驱动 力 分 析 ， 
也 多 基于 区 域 整 体 性 ,而 缺乏 区 域内 部 之 间 的 差异 
性 对 比 ,特别 是 对 于 流域 生态 环境 质量 的 预测 研究 
鲜 见 报道 。 

鉴于 此 ,本 文 基于 Landsat 遥感 影像 .DEM 数 
据 、 气 象 数 据 和 人 口 经 济 数据 等 ,采用 RSEI 评 价 指 
数 , 对 延 河 流域 1998 一 2016 年 的 生态 环境 质量 进行 
时 空 测评 ,对 生态 环境 质量 下 降 的 区 域 进行 驱动 力 
探究 ,并 尝试 利用 CA-Markov 模 型 ,对 延 河流 域 2022 
年 和 2028 年 的 生态 环境 发 展 格局 进行 模拟 和 预测 ， 
为 流域 生态 安全 规划 重点 和 环境 治理 保护 的 制定 
提供 科学 依据 。 
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1 研究 区 概况 与 方法 


1.1 研究 区 概况 

延 河 流域 (108°45’~110°28'E,36°23’~37°17’'N) 
位 于 我 国 西北 内 陆 半 湿润 向 半 干 旱 过 渡 地 区 (图 
1) ,属于 我 国 黄土 夺 陵 沟 密 区 第 二 副 区 。 延 河 是 黄 
河 的 一 级 支流 ,干流 河 长 约 为 286.9 km, 发源 于 榆林 
市 靖 边 县 天 赐 湾 乡 的 周 山 , 流 经 延安 市 的 志 丹 县 、 
安 塞 县 FAK MERE RAILCAR, yEy 
流域 地 势 西 北 高 东南 低 , 属 暧 温带 向 中 温带 过 渡 的 
大 陆 性 季风 气候 区 ,多 年 平均 气温 为 9.2 % ,降水 量 
为 500 mm ,土壤 类 型 主要 为 黄 绢 土 。 流 域 水 土 流失 
较为 严重 ,多 年 的 输 沙 量 高 达 0.83x108t ,侵蚀 模 数 
ZJ 1.12108 t+ km 
1.2 数据 来 源 及 预 处 理 

人 研究 选用 云 量 低 于 10% 且 季节 相同 的 Landsat 
影像 数据 , 辅 以 DEM .降雨 .温度 .人口 和 GDP 数 
据 。 各 数据 人 处理 的 具体 方法 如 表 1 所 示 。 其 中 , 辐 
射 定 标 是 为 了 将 遥感 影像 的 DN 值 转换 为 传感器 处 
的 反射 率 , 而 大 气 校正 则 为 了 消除 大 气 和 光照 等 因 
素 对 地 物 反 射 的 影响 。 
1.3 水 体 掩 膜 获 取 

获取 遥感 生态 指数 需要 使 用 主 成 分 变换 方法 ， 
将 代表 环境 状态 的 4 个 指标 干 度 (NDBSI) 湿度 
(Wet) . 绿 度 (NDVI) 和 热度 (LST) HET AA. FA 
于 大 面积 水 域 会 对 真实 地 面 湿 度 以 及 主 成 分 的 载 


105°E 


108°E 111°E 


(a) 延 河流 域 位 置 J 


39°N 


36°N 


33° N 


37°0'N 


36°30'N 


荷 产 生 影响 ,因此 需 先 剔除 大 面积 水 域 再 进行 主 成 
分 变换 ”。 
研究 选用 改进 的 归 一 化 差异 水 体 指数 (MND- 
WI) 获 取 流 域内 的 水 体 掩 膜 。 该 指数 是 将 归 一 化 差 
异 水 体 指数 公式 中 的 近 红 外 波段 改 为 了 中 红外 波 
段 ,以 此 来 提高 水 体 和 建筑 物 的 分 类 精度 1。 
MNDWI = (Green - MIR)/(Green + MIR) (1) 
式 中 :Green 为 绿 波 段 反 射 率 ;MIR 为 中 红外 波段 反 
射 率 。 分 别 对 应 LandsatS TM 的 第 2. 5 波段 反射 率 
All Landsat8 OLI 的 第 3.6 波 段 反 射 率 。 
1.4 RSEI 计 算 
RSEI 选 用 主 成 分 分 析 方 法 ,将 绿 度 .湿度 FRE 
和 热度 4 个 指标 的 信息 集中 到 第 一 主 成 分 上 ,其 是 
根据 各 指标 在 主 成 分 上 的 贡献 度 和 自身 的 性 质 来 
确定 权重 值 ,减少 了 人 为 赋值 的 误差 。 为 了 避免 水 
体 对 主 成 分 载荷 的 影响 和 指标 量 纲 不 统一 对 权重 
的 影响 ,在 进行 主 成 分 分 析 之 前 需 对 各 指标 进行 如 
下 两 步 处 理 :(1) 利用 水 体 掩 膜 去 除 水 体 ,(2) 对 4 
个 指标 进行 归 一 化 。 
指标 归 一 化 公式 如 下 : 
NL,=(1,—7,)/(T ~ Lin) (2) 
式 中 :NL 为 各 指标 的 归 一 化 结果 ;1 为 各 指标 在 第 i 
像 元 的 值 ;Li 和 ,分 别 为 各 指标 的 最 大 值 和 最 
小 值 。 
RSEI 计 算 公 式 如 下 : 
RSEI, = 1 - PC1[ f (NDVI, Wet, LST, NDBSI)] (3) 
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图 1 延 河流 域 地 理 位置 
Fig. 1 Geographical location of Yanhe River Basin 
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表 1 研究 数据 和 处 理 方法 


Tab. 1 Research data and processing methods 


数据 数据 来 源 时 间 空间 分 辨 率 /m 处 理 方法 
影像 Landsat TM 地 理 空间 数据 云 1998 年 6 月 16 日 30 辐射 定 标 、 大 气 校正 影像 拼 接 和 裁剪 等 预 处 理 
数据 http://www.gscloud.cn/ ”2004 年 6 月 14 日 
2010 年 6 月 30 日 
Landsat OLI 地 理 空间 数据 云 2016 年 6 月 17 日 30 辐射 定 标 、 大 气 校正 .影像 拼接 和 裁剪 等 预 处 理 
http://www.gscloud.cn/ 
辅助 DEM 地 理 空间 数据 云 2009 年 30 投影 变换 .拼接 和 裁剪 
数据 http://www.gscloud.cn/ 
逐日 降雨 和 中 国 科学 院 资源 环境 1996 年 .1998 年 30 选取 研究 区 内 、 外 共 27 个 站 点 数据 ,计算 各 站 点 
温度 数据 科学 与 数据 中 心 2000 年 .2002 年 年 均 温度 和 年 降水 量 , 采 用 反 距离 插值 法 将 计算 
hitps://www.resde.cn/ 2004 4F „2006 4E 得 到 的 数据 搬 值 为 空间 分 辩 率 与 遥感 影像 一 至 
2008 年 .2010 年 的 栅 格 数据 , 求 取 均 值 
2012 年 .2014 年 
2016 年 .2018 年 
人 口 和 GDP 陕西 省 统计 年 鉴 1998 一 2019 年 - 对 延 河流 域 各 县 人 口 数据 和 GDP 数 据 进行 统计 


注 : 表 中 -表示 无 相应 属性 。 


RSEI= (RSEI, ~ RSEI,,,)/(RSEI,,,, -RSEI,,,) (4) 
式 中 :PC1 为 4 个 归 一 化 后 的 指标 通过 主 成 分 分 析 
方法 所 得 到 的 第 一 主 成 分 结果 ,为 了 方便 分 析 ,对 
其 进行 正 负 值 转 置 ” ,RSEI 为 转 置 后 结果 ,4 个 指标 
分 别 参 考 徐 涵 秋 "等 的 研究 成 果 , 具 体 计算 公式 
如 表 2 所 示 。RSEIu. 和 RSEIus 分 别 为 RSEl 的 最 大 
值 和 最 小 值 。 

1.5 CA-Markovy 模 型 
CA-Markov 是 一 个 可 以 基于 像 元 进行 时 空 演变 


模拟 预测 的 模型 , 它 集成 了 元 胞 自动 机 (CA ) 模 型 的 
空间 模拟 优势 和 马尔 科 夫 (Markov) 模 型 的 时 间 预 测 
优势 。 此 模型 普遍 适用 于 土地 利用 变化 的 模拟 和 
预测 , 因 生 态 环 境 质 量 也 是 具有 高 度 空 间 自 相关 性 
的 栅 格 数据 ,因此 人 研究 尝试 使 用 CA-Markov 模 型 对 
延 河流 域 的 生态 环境 质量 进行 预测 。 

CA 模型 是 一 种 时 间 、 空 间 和 状态 均 离 散 的 局 部 
网 格 动 力学 模型 ,该 模型 具有 模拟 复杂 系统 时 空 演 
化 过 程 的 能 力 “。 其 公式 如 下 : 


表 2 各 指标 公式 


Tab.2 Calculation formula of each index 


HER 公式 参数 说 明 
NDVI NDVI=(NIR-Red)/(NIR+Red) Red 为 红 波段 反射 率 ;NIR 为 近 红 外 波段 反射 率 。 分 别 可 对 应 
Landsat5 TM 的 3. 4 波段 反射 率 和 Landsat8 0LI 的 4.5 波 段 反 射 率 
Wet Wet(TM) = 0.0315p, + 0.202 Ip, + 0.3102p, + p (1.2.3.4.5.0) 分 别 代表 Landsat5 TM 的 1.2.3.4.5.7 波 段 反 射 率 
0.1594p, - 0.6806p, — 0.6109, All Landsat8 0LI 的 2.3、4、5、6.7 波 段 反射 率 
Wet(OLI)=0.1511p, + 0.1973p, + 0.3283p, + 
0.3407p, - 0.7117p, — 0.4559p。 
NDBSI NDBSIe (SI + IBI) p: (1.2.3.4.3) 分 别 代表 Landsat5 TM 波 段 1.2.3.4.5 的 反射 率 和 
2 Landsat8 OLI 波 段 2.3.4.5.6 的 反射 率 
SI= [Cp; +p;) 一 (ps tp)\/[(es +p;) + (ps +p,)| 
pie! 2 -| Pa , P | 
PstPs [PstP, Ps+p; 
2p; + ps + P2 | 
pstps |Pp3tps  PrtPs 
LST DN 为 Landsat 数 据 像 元 的 灰 度 值 ;gain 和 bias 分 别 为 波段 增益 值 和 


rsT=/7,l + (7)/p)e|-237.15 


< K, 
T,= Kj inf + 1) 


L, = gain X DN + bias 


A EEIE ; Lof RÆ TM/TIRS ATAN BS HIETE ; K PIK AEREA 
7 为 亮度 温度 ;A 为 TM 影像 第 6 波段 的 中 心 波 长 11.5 um Fil OLI AK 
据 的 第 10 波 段 中 心 波长 10.9 pm;p 等 于 1.438x10”mK;s 为 比 辐 
射 率 


注 :NDVI 代 表 绿 度 指标 ;Wet 代 表 湿 度 指标 ;NDBSI 代 表 干 度 指 标 ;LST 代 表 热 度 指标 。 
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Seayn) (5) 
式 中 :5 为 元 胞 状态 ;f 为 局 部 空间 的 元 胞 转化 规则 ; 
titl 为 各 元 胞 时 刻 ;n 为 2 个 时 刻 的 间隔 时 间 。 

Markov 模 型 认为 事物 的 变化 服从 马尔 可 夫 过 

程 ,其 认为 事物 的 后 一 个 状态 仅 与 前 一 时 刻 状态 和 
转移 概率 矩阵 有 关 , 其 公式 如 下 : 

Xa = Xo XP (6) 
式 中 ;Xow 分 别 为 :+1\t 时 刻 的 状态 ;Pj; 为 转移 概 
率 和 矩阵 。 


2 结果 与 分 析 


2.1 各 指标 主 成 分 分 析 

由 表 3 可 知 ,4 期 的 第 一 主 成 分 特征 值 贡献 率 均 
在 80 狗 左右 , 且 各 指标 的 贡献 率 均 仅 在 第 一 主 成 分 
上 有 着 相同 的 正 负 分 布 ,具体 为 代表 植被 覆盖 度 的 
NDVI 和 代表 生态 环境 湿度 的 Wet 在 第 一 主 成 分 上 
的 贡献 率 为 正 值 ,代表 地 表 温 度 的 LST 和 代表 建筑 


和 裸 土 所 用 地 硬化 程度 的 NDBSI 为 负 值 ,表明 ND- 
VI 和 Wet 对 生态 环境 的 改善 起 促进 作用 ,LST 和 
NDBSI 对 生态 环境 的 改善 起 抑制 作用 2 ,这 与 现实 
中 4 个 指标 对 生态 环境 的 影响 相符 合 。 因 此 ,将 选 
用 第 一 主 成 分 的 结果 进行 RSEI 计 算 。 
2.2 流域 生态 环境 质量 时 空 变 化 

在 时 间 尺 度 上 ,对 各 年 份 4 个 指标 和 RESI 进 行 
统计 ( 表 4) ,由 表 4 可 以 看 出 :(1) RSEI 指 数 均值 不 
断 增加 ,从 1998 年 的 0.392 增 长 为 2016 年 的 0.530， 
表明 流域 生态 环境 质量 整体 有 所 改善 ,特别 是 
2010 一 2016 年 期 间 ,RSEI 指 数 均值 增长 了 20.24%， 
生态 环境 质量 提升 较 快 。(2) 1998 一 2016 年 , NDVI 
均值 呈 增 长 趋势 ,与 韩 右 等 的 研究 成 果 相 符 。 
Wet 均 值 整体 呈 上 升 趋势 。LST 均 值 仅 在 2016 年 有 
所 降低 。NDBSI 呈 下 降 趋势 ,可 能 是 由 于 近 几 年 延 
安 市 实行 了 退耕 还 林政 策 , 大 面积 耕地 转化 为 林地 
使 得 裸 土 指数 有 所 降低 -(3) 对 比 各 指数 在 PC1 上 
荷载 的 绝对 值 , 发 现 NDVI、Wet 和 NDBSI 近似 相等 


表 3 1998 年 .2004 年 .2010 年 和 2016 年 各 指标 主 成 分 分 析 结 
Tab. 3 Principal component analysis results of each index in 1998, 2004, 2010 and 2016 


年 份 参数 第 一 主 成 分 第 二 主 成 分 第 三 主 成 分 第 四 主 成 分 
1998 NDVI 0.521 0.420 -0.383 0.636 
Wet 0.525 -0.102 0.834 0.139 
NDBSI -0.530 -0.340 0.166 0.759 
LST -0.414 0.835 0.362 0.006 
特征 值 0.137 0.022 0.008 0.000 
特征 值 贡 献 率 /9 81.640 13.350 4.850 0.160 
2004 NDVI 0.493 0.418 -0.478 -0.595 
Wet 0.514 0.119 0.834 -0.161 
NDBSI -0.534 -0.242 0.213 -0.782 
LST -0.456 0.868 0.175 0.091 
特征 值 0.126 0.021 0.009 0.001 
特征 值 贡献 率 /% 79.950 13.150 5.680 0.870 
2010 NDVI 0.548 0.328 -0.388 0.665 
Wet 0.547 0.003 0.837 0.036 
NDBSI -0.541 -0.223 0.322 0.744 
LST -0.329 0.918 0.214 -0.057 
特征 值 0.138 0.019 0.006 0.001 
特征 值 贡献 率 /% 84.050 11.610 3.570 0.770 
2016 NDVI 0.565 0.298 -0.548 0.541 
Wet 0.502 0.061 0.819 0.273 
NDBSI -0.558 -0.216 0.093 0.796 
LST -0.344 0.928 0.144 -0.006 
特征 值 0.167 0.021 0.004 0.001 
特征 值 贡献 率 /% 86.720 10.710 2.280 0.290 
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表 4 各 年 份 4 个 指标 和 RSEI 的 统计 值 


Tab.4 Statistical values of 4 indicators and RSEI for each year 


年 份 参数 NDVI Wet NDBSI LST RSEI 
1998 均值 0.322 0.478 0.685 0.507 0.392 
对 PC1 和 荷载 值 0.521 0.525 -0.530 -0.414 一 
2004 均值 0.323 0.507 0.675 0.513 0.410 
对 PC1 和 荷载 值 0.493 0.514 -0.534 -0.456 = 
2010 均值 0.471 0.503 0.619 0.539 0.441 
对 PC1 和 荷载 值 0.548 0.547 -0.541 -0.329 一 
2016 均值 0.581 0.565 0.504 0.420 0.530 
对 PC1 荷载 值 0.565 0.502 -0.558 -0.344 - 


注 : 表 中 -表示 无 相应 数值 。 


日 大 于 LST, 表 明 该 流域 生态 环境 易 受 NDVI、Wet 和 
NDBSI 的 影响 。 

参考 《生态 评价 技术 规范 ) 中 的 生态 环境 分 级 
标准 和 以 往 学 者 的 研究 成 果 ” ,将 计算 得 到 的 RSEI 
值 以 0.2 为 间隔 划分 为 很 差 较 差 .中 等 .良好 和 优 
恨 5 个 等 级 (图 2) ,并 对 各 等 级 面积 进行 统计 ( 表 
5), 由 表 5 可 以 看 出 :(1) 空间 尺度 上 ,流域 生态 环 
境 质量 基础 较 差 ,总体 分 布 特点 表现 为 南部 优 于 北 
部 ,主要 原因 是 南部 水 热 条 件 组 合 较 好 。 生 态 环 境 


(a) 1998 年 


Ta. TA, 


(€) 20104 


质量 以 较 差 和 中 等 等 级 为 主 ,主要 分 布 于 河流 两 旁 
的 山谷 地 带 ,平均 面积 占 比 为 62.4%。 优 良 等 级 面 
积 占 比 最 少 , 主要 分 布 于 流域 南部 的 山区 地 带 , 平 
均 面 积 占 比 仅 为 8.1%。(2) 生态 环境 质量 从 东南 向 
西北 逐渐 改善 。 较 差 和 良好 等 级 的 面积 占 比 变化 
显著 ,其 中 较 差 等 级 面积 占 比 减少 了 19.3%, 良 好 等 
级 面积 占 比 增加 了 16.5%。(3) 2012 年 4 月 流域 开展 
的 “ 治 沟 造 地 ”工程 中 ,对 生态 环境 质量 影响 较 大 ， 
图 2 中 显示 ,到 2016 年 流域 中 部 新 区 的 生态 环境 质 


(b) 2004 年 


(d) 2016 年 


a, A. 


RSEI 


Mkt Re DS May MR 


0 15km 
上 一 


图 2 延 河 流域 1998 年 .2004 年 .2010 年 和 2016 年 遥感 生态 指数 
Fig. 2 Remote Sensing ecological index of Yanhe River Basin in 1998, 2004, 2010 and 2016 
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表 5 流域 1998 一 2016 年 遥 


态 指数 各 等 级 面积 统计 


Tab.5 Area statistics of each grade of remote sensing pa index in watershed from 1998 to 2016 


RSEI 1998 4 2004 4E 20104 2016 4 均值 

等 级 ”面积 km ” 占 比 /% 面积 km 。 占 比 /%6 面积 km He% 面积 km 占 比 /% 面积 km 占 比 /% 
很 差 1240.430 16.2 1344.149 17.6 948.261 12.4 687.401 9.0 1055.060 13.8 
较 差 3123.612 40.9 2731.561 35.7 2563.76 33.5 1650.651 21.6 2517.396 32.9 
中 等 。 2177.107 28.5 2167.987 28.3 2428.126 31.8 2229.879 29.2 2250.775 29.5 
良好 706.533 9.3 878.767 11.5 1252.57 16.4 1969.7 25.8 1201.892 15.7 
优良 390.561 5.1 528.387 6.9 452.606 5.9 1103.360 14.4 618.729 8.1 
合计 7638.244 100.00 7650.851 100.00 7645.323 100.00 7640.991 100.00 7643.852 100.00 

量 等 级 变 为 很 差 。 和 志 丹 县 等 区 域 生 态 环境 质量 恶化 ,其 余 各 县 的 


23 延 河流 域 生态 环境 质量 动态 监测 


为 了 定量 研究 4 期 延 河流 域 生态 环 境 质 量变 化 
的 时 空 分 布 特征 和 变化 趋势 ,采用 差 值 法 对 4 期 的 


RSEI 进行 


空间 对 比 ,并 划分 为 变 好 、 不 变 和 变 差 3 
个 类 别 (图 3)。 


由 表 6 可 以 看 出 ,流域 生态 环境 质量 


保持 稳定 的 面积 占 比较 大 ,同时 有 所 改善 。 而 且 差 


值 的 等 级 数 为 1 的 面积 最 大 ,说 明 流 域 生态 环境 质 
量 总 体 表现 为 循序 渐进 地 向 高 一 等 级 转变 
由 图 3 可 以 看 出 ， ge 、 宝 塔 


(a) 1998 一 2004 年 


(c) 2010 一 2016 年 


T Ra, 


Fig. 3 Spatiotemporal distribution of the change of ecological environmental quality in Yanhe River Basin 


eR eae ete ;主要 因 前 者 多 
位 于 城市 周边 ,建设 用 地 扩张 一 定 程度 上 影响 了 生 
态 环境 质量 ,后 者 生态 退耕 政策 的 实施 起 到 了 积极 
的 作用 。2004 一 2010 年 ,流域 西部 和 中 部 区 域 生态 
环境 质量 有 所 改善 ,西北 部 以 及 南部 区 域 生 态 环境 

质量 变 差 ,可 能 与 这 个 时 期 流域 土壤 侵蚀 强度 增加 


AR", 
面 改善 


(b) 2004—20104F 


T Wa, 


(d) 1998 一 2016 年 


O ELA mee DAS mee 0 1 


3 延 河流 域 


E 态 环境 质量 变化 时 空 


分 布 


2010 一 2016 年 ,流域 生态 环境 质量 得 到 全 
。 这 也 表明 ,政策 因素 对 于 生态 环境 的 质量 
具有 重要 的 驱动 作用 。 
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表 6 延 河 流域 生态 等 级 变化 统计 


Tab.6 Statistics on the change of ecological grades in Yanhe River Basin 


类 别 关 什 1998 一 2004 年 2004 一 2010 年 2010—2016 年 1998 一 2016 年 
面积 /km 占 比 /% 面积 /km? 占 比 /% 面积 /km 占 比 /% 面积 /km 占 比 /% 
变 差 -4 0.118 0.000 0.572 0.010 2.343 0.030 1.577 0.020 
-3 2.410 0.030 6.343 0.080 12.629 0.170 13.319 0.170 
-2 68.338 0.890 67.848 0.890 64.073 0.840 116.252 1.520 
=1 1198.591 15.690 1224.513 16.020 506.432 6.630 891.650 11.690 
不 变 0 4499.537 58.910 4238.542 55.460 3571.693 46.770 2596.093 34.020 
变 好 1 1764.221 23.100 1737.903 22.740 2894.729 37.910 2328.290 30.510 
2 99.059 1.300 286.970 3.750 518.229 6.790 1085.440 14.230 
3 5.306 0.070 69.857 0.910 63.083 0.830 449.429 5.890 
4 0.428 0.010 10.327 0.140 3.549 0.050 148.193 1.940 


2.4 驱动 力 因 子 分 析 
241 影响 生态 环境 质量 交差 的 因素 分 析 RSEI 
的 4 个 指标 分 别 代 表 着 流域 植被 覆盖 、 土 地 退化 、 土 
二 湿度 和 地 表 温 度 情况 ,是 指 征 生态 环境 质量 的 重 
要 指标 。 因 此 ,针对 生态 环境 质量 恶化 区 域 , 分 析 
各 指标 变化 对 RSEI 变 化 的 驱动 力 , 有 助 于 更 加 清晰 
地 了 解 生态 环境 质量 变 差 的 原因 ,以 便 为 流域 生态 
环境 改善 提供 科学 依据 。 

使 用 GIS 的 空间 分 析 方 法 ,可 视 化 表达 各 指标 


c 区 域 主要 受 Wet 和 NDBSI 变 化 的 影响 ,d 区域 主要 
受 Wet、NDBSI 和 LST 变 化 的 影响 ,而 e 区 域 主要 受 
NDVI 和 Wet 变 化 的 影响 。 影 响 土壤 湿度 的 因素 有 
很 多 ,但 起 重要 作用 的 为 降水 量 ag Ac at Bl HES 
型 ? ,土地 退化 则 与 气候 变化 和 土壤 侵蚀 有 着 密切 
的 关系 '。 分 析 流 域 多 年 的 气候 变化 情况 (图 5)， 
可 以 发 现 气 温和 降水 量 多 年 变化 幅度 较 大 ,这 在 一 
定 程 度 上 会 加 剧 流域 的 水 土 流失 问题 ,从 而 导致 生 
态 环 境 质量 变 差 。 因 此 ,a、b.c 和 d 区 域 生态 环境 质 


1998 一 2016 年 的 变化 情况 (图 4) ,划分 生态 环境 质 
量变 差 区 域 , 由 图 4 可 以 看 出 ,其 主要 分 布 于 a、be、 
d 和 e 区 域 ,并 采用 SPSS 工具 建立 各 区 域 的 多 元 线 
性 回归 模型 ( 表 7), 由 表 7 可 以 看 出 ,各 区 域 对 RSEI 
等 级 降低 起 重要 作用 的 指标 存在 差异 。 其 中 ab 和 


(a) NDVI 


(b) RSEI 


图 例 
m 变 差 
口 不 变 
mm 变 好 


(d) LST 


图 例 
四 变 差 
口 不 变 
me 变 好 


量变 差 较 多 受 自 然 因素 影响 。 植 被 覆盖 度 的 变化 
与 人 为 因素 和 自然 因素 都 有 着 密切 的 关系 ,而 延 河 
流域 人 地 矛盾 问题 突出 ,人 口 数 量 和 经 济 收入 多 年 
呈 增 长 趋势 (图 5) ,加 速 了 流域 城市 的 开发 ,并 且 e 
区 域 为 延安 市 新 区 ,因此 该 区 域 生态 质量 变 差 较 多 


(c) Wet 


图 例 


me 变 差 
口 不 变 
m 变 好 


(e) NDBSI 


图 例 
m 变 差 
口 不 变 
mm 变 好 


30 km 
a | 


图 4 1998 一 2016 年 延 河 流域 各 指标 变化 
Fig. 4 Changes of indicators in Yanhe River Basin from 1998 to 2016 
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表 7 多 元 线性 回归 方程 


Tab.7 Multiple linear regression equation table 


区 域 样本 点 数 多 元 线性 回归 方程 R R P 
a 158 RSEL=1.688+0.147NDVI+0.347W et—0.363NDBSI-0.009LST; 0.789 0.622 <0.01 
b 300 RSEL=2.332+0.086NDVI;+0.199W et—0.238NDBSI-0.220LST; 0.585 0.342 <0.01 

250 RSEL=2.520+0.012NDVI;+0.195 Wet—0.384NDBSI-0.032LST; 0.643 0.413 <0.01 
d 534 RSEL=1.335+0.150NDVI;+0.368 Wet—0.188NDBSI-0.090LST; 0.641 0.411 <0.01 
e 100 RSEL=1.201+0.385NDVI+0.263 Wet-0. 13 1NDBSI-0.112LST; 0.772 0.596 <0.01 


iE: RSEL,NDVL Wet: NDBSL LST 为 1998 一 2016 年 各 指标 值 差 划分 等 级 后 的 结果 值 (其 划分 方法 与 图 3 相同 ) , 取 值 为 1.2.3; 其 中 ,等 级 变 差 
取 值 为 1; 等 级 不 变 取 值 为 2; 等 级 变 好 取 值 为 3; 下 标 i 为 等 级 数 , 取 值 为 1,2,3。 


(a) 年 降水 量 和 年 平均 温度 

— 年 降水 量 /mm 一 年 平均 温度 /'C 
700 
600 
500 


年 降水 量 /mm 
年 平均 气温 /'C 


0 1 1 1 1 1 
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(b) 人 口 和 GDP 
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图 5 延 河流 域 多 年 气候 和 社会 经 济 情况 


Fig. 5 Climatic and socio-economic conditions in Yanhe River Basin 


受 人 为 因素 影响 。 

2.4.2 DEM 与 生态 环境 质量 的 关系 ， 基 于 GIS 空 间 
统计 功能 分 析 DEM 与 RSEI 等 级 之 间 的 关系 ,发 现 
海拔 低 于 600mm 的 区 域 ,生态 环境 质量 为 很 差 . 较 差 
和 中 等 的 等 级 面积 占 比 超过 90% ,优良 等 级 占 比 几 
乎 为 0( 图 6)。 在 海拔 600~1800 m, 随 着 海拔 的 升 
高 ,流域 很 差 等 级 面积 占 比 增加 ,海拔 1600~1800 m 


E 很 差 
1600~1800 


eRe m 中 等 


很 差 和 较 差 等 级 面积 占 比 超过 了 80%。 参 考 《 生 态 
环境 状况 评价 技术 规范 》 中 对 各 等 级 的 描述 ,以 具 
备 人 类 居住 条 件 的 中 等 .良好 和 优良 等 级 为 划分 标 
准 , 发 现 生态 环境 质量 以 1200 m 为 分 界 点 ,473~ 
1200 m 随 着 海拔 的 升 高 具备 人 类 居住 条 件 的 面积 
占 比 增加 ,1200~1800 m 随 着 海拔 的 升 高 面积 占 比 
减少 ,表明 延 河 流域 在 较 高 和 较 低 海拔 区 域 的 生态 


eRe 加 优良 


1400~1600 M 


REE 


xe 800~1000 


一 
z 1000-1200 fl 


6oo-soo Jill 
473-600 ii 和 有 国 
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各 等 级 面积 占 比 /% 


图 6 不 同 海拔 不 同等 级 生态 质量 所 占 比 
Fig.6 Proportions of different levels of ecological quality at different altitudes 
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(a) 2022 年 生态 环境 质量 模拟 值 
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(b) 2028 年 生态 环境 质量 预测 值 


图 7 2022 年 (a) 模 拟 与 2028 年 (b) 预 测 结 


Fig.7 Simulated results in 2022 (a) and predicted results in 2028 (b) 


环境 质量 较 差 ,在 海拔 800~1200 m 区 域 生 态 环境 质 
量 较 好 ,大 面积 适宜 人 类 居住 。 
2.5 CA-Markov 模 型 预测 

基于 1998 年 .2004 年 和 2010 年 的 RSEI 等 级 计 
算 结 果 ,运用 CA-Markov 模型 ,模拟 2010 年 和 2016 
年 的 RSEI 等 级 ,并 评价 模型 的 适用 性 。 经 过 计算 ， 
2010 年 和 2016 年 的 模拟 结果 与 计算 结果 之 间 的 
Kappa 系数 分 别 为 0.6863 和 0.6102 ,属于 高 度 一 致 ， 
表明 该 模型 在 延 河 流域 生态 环境 质量 的 预测 中 具 
有 一 定 的 适用 性 。 因 此 ,应 用 此 模型 对 流域 2022 年 
的 生态 环境 情况 进行 模拟 ,并 基于 模拟 结果 对 2028 
年 进行 预测 (图 7)。 

经 统计 ,2022 年 等 级 很 差 较 差 中等、 良好 和 
优良 面积 各 占 比 为 6%、15%、24%、30% 和 25% ,2028 
年 为 5% 、11%、18%、31% 和 35%。 相 比 于 2016 年 ， 
很 差 较 差 和 中 等 等 级 面积 减少 ,其 余 各 等 级 面积 
均 增 加 ,其 中 等 级 较 差 ,优良 面积 变化 最 大 。 到 
2028 年 延 河 流域 生态 环境 的 改善 会 持续 向 西北 方 
向 延伸 ,但 由 于 流域 西北 部 地 人 处 半 干 旱地 区 ,降水 
稀少 、 植 被 覆盖 度 较 低 .生态 环境 本 底 值 较 差 。 
此 ,未 来 对 流域 生态 环境 治理 的 重点 仍 应 为 西北 部 
地 区 。 


3 讨论 


目前 ,遥感 生态 指数 主要 应 用 于 对 城市 和 矿区 
生态 环境 质量 的 评价 ,对 流域 生态 环境 的 监测 ,多 
集中 于 半 湿 润 区 和 湿润 地 区 ,在 半 湿 润 区 向 半 干 旱 
区 过 度 地 带 的 应 用 较 少 。 因 此 ,应 用 RSEI 对 生态 环 


境 本 底 脆弱 、 人 地 关系 矛盾 较为 突出 的 延 河 流域 进 
行 生态 环境 监测 评价 具有 一 定 挑战 。 人 研究 表明 ， 
RSEI 值 估算 结果 与 延 河流 域 环境 治理 工程 实施 前 
后 ,流域 生态 环境 质量 的 改善 状况 相 吻 合 。 同 时 ， 
一 些 学 者 的 研究 结果 表明 ,延安 市 植被 覆盖 度 增加 
的 总 体 迁移 趋势 是 向 北向 西 , 北 部 植被 覆盖 度 仍 比 
较 低 ” ,这 也 与 本 人 研究 所 监测 的 结果 一 致 ,表明 
RSEI 对 延 河 流域 具有 较 好 的 适用 性 。 

基于 此 ,研究 利用 RSEI 对 延 河 流域 近 20 a 的 生 
态 环 境 进 行 遥 感 监 测 ,结果 表明 延 河流 域 的 生态 环 
境 本 底 较 弱 ,但 发 展 较 好 。 海 拔 较 高 地 区 的 生态 环 
境 质 量变 差 主 要 是 受 自然 因素 影响 ,其 中 土壤 干 湿 
度 和 地 表 温 度 的 影响 较为 突出 ,而 海拔 适宜 地 区 的 
生态 环境 质量 变 差 主要 是 受 人 类 活动 影响 ,例如 城 
市 建设 用 地 发 展 等 。 柯 锐 鹏 等 ”研究 发 现 , 当 城 镇 
化 发 展 进 入 高 级 阶段 ,城镇 化 可 以 保持 甚至 改善 生 
态 环境 质量 ,这 就 需要 人 们 对 城市 进行 合理 规划 。 
因此 ,需要 注重 流域 西北 部 地 区 生态 环境 的 综合 治 
理 以 及 中 部 地 区 城市 规模 的 适度 发 展 。 

文中 根据 各 个 指标 的 变化 情况 ,分析 了 流域 生 
态 环境 质量 变 差 区 域 的 驱动 因子 ,但 研究 中 着 重 考 
虑 了 各 分 区 的 主导 驱动 因素 ,缺乏 对 各 因子 之 间 的 
交互 作用 的 综合 考虑 ,也 未 实现 各 驱动 因子 的 定量 
人 研究。 因此 ,后 续 的 研究 中 , 拟 采 用 地 理 探测 带 方 
法 ,对 流域 的 生态 环境 变化 的 驱动 因子 进行 定量 分 
析 。 同 时 ,在 对 CA-Markov 模 型 的 可 适用 性 进行 检 
验 时 ,2 期 模拟 结果 的 Kappa 系数 虽 符 合 高 度 一 致 
性 ,但 也 存在 一 定 误差 ,这 与 影像 本 身 的 不 确定 性 
以 及 因素 选取 的 全 面 性 都 有 关系 ,也 期 望 后 续 能 得 
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到 改进 。 
4 结论 


利用 Landsat 影 像 数 据 , 辅 以 DEM 气象 数据 和 
人 口 经 济 数据 ， me 2004 年 .2010 


年 和 2016 年 的 生态 环 行 分 析 和 评价 ， tie 
bi CA-Markov 模型 aioe 2022 es 2028 年 的 生态 环 
进行 预测 。 得 到 如 下 结论 : 


(1) es E 质量 评价 中 
具有 较 好 的 适用 性 ， ee P 
环境 质量 监测 中 可 进行 推广 应 用 。 
(2) 延 河流 域 生 态 基 础 较 差 ,但 近 20 a 来 生态 
环境 质量 在 不 断 提 升 。RSEI 均 值 从 1998 年 的 0.392 
增长 到 2016 年 的 0.530, 且 空间 上 生态 环境 质量 从 
东南 向 西北 逐渐 改善 。 
(3) 自然 因素 和 人 为 因素 对 流域 生态 环境 变化 
的 影响 存在 空间 差异 性 。 其 中 海拔 1200~1800 m 的 
区 域 , 生 态 环境 质量 变 差 主 要 受降 雨量 .气温 变化 
等 自然 因素 的 影响 ; 低 于 1200 m 的 区 域 ,主要 受 城 


市 扩张 等 人 类 活动 的 影响 。 在 海拔 800~1200 m, 受 
自然 因素 、 人 类 活动 的 影响 均 较 少 ,生态 环境 质量 


普遍 高 于 其 他 区 域 。 因 此 ,高 海拔 地 区 土地 资源 的 
合理 优化 以 及 中 低 海拔 地 区 城市 化 进程 的 适度 发 
展 是 值得 关注 的 问题 

(4) CA-Markov 模型 的 预测 结果 表明 ,2028 年 
延 河 流域 生态 环境 质量 会 得 到 大 面积 改善 ,良好 和 
优良 等 级 面积 占 比 总 和 可 达到 66% ,但 西北 部 还 是 
以 很 差 和 较 差 等 级 为 主 。 因 此 ,未 来 对 于 生态 本 底 
脆弱 的 西北 部 地 区 仍 需 加 强 保护 和 治理 。 
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Monitoring and prediction of ecological environmental quality in the Yanhe 
River Basin based on the remote sensing ecological index 


WANG Lixia, ZHAO Rui, LIU Zhao’, ZHANG Shuangcheng', 
KONG Jinling', YANG Yun' 


(1. School of Geological Engineering and Geomatics, Chang’ an University, Xi an 710054, Shaanxi, China; 


2. School of Water Conservancy and Environment, Chang’ an University, Xi an 710054, Shaanxi, China) 


Abstract: Based on Landsat remote sensing images of the Yanhe River Basin in 1998, 2004, 2010, and 2016, the 
remote sensing ecological index (RSEI) was calculated to analyze and evaluate temporal and spatial changes of 
environmental quality in the Yanhe River Basin. The CA-Markov model was also used to simulate and predict the 
ecological environment in 2022 and 2028. Results show that (1) the RSEI has good applicability in the Yanhe 
River Basin and can be used to monitor and evaluate temporal and spatial characteristics of the ecological 
environment in the basin. (2) For the time scale, the mean RSEI value changed from 0.392 in 1998 to 0.530 in 
2016 and showed an increasing trend year by year. On a spatial scale, the RSEI gradually increased from 
southeast to northwest. RSEI grades were divided by 0.2 intervals and the total area of RSEI grades was mainly 
poor and moderate. These grades accounted for a total of 62.4%, which was distributed mainly in the valley of the 
basin, while the high grade area accounted for the least, at approximately 8.1%, which was distributed mainly in 
the mountainous region of the southern part of the basin. (3) There is a large spatial difference between the effects 
of natural and human factors on changes in environmental quality in the watershed. In the area with an altitude of 
1200-1800 m, the change of eco-environmental quality was affected mainly by natural factors. However, areas 
below 1200 m are affected mainly by human activities. (4) It is predicted that by 2028, the ecological 
environment quality of the basin will continue to improve, and the area of good and high grade will account for 
66%, but there are still big problems in the northwest of the basin because of the fragile ecological environment 
background. 

Keywords: Yanhe River Basin; remote sensing ecological index; ecological environment quality; principal com- 
ponent analysis; CA-Markov model 


